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Abstract 

Background: Some species of water frogs originated from hybridization between different species. Such hybrid 
populations have a particular reproduction system called hybridogenesis. In this paper we consider the two species 
Pelophylax ridibundus and Pelophylax lessonae, and their hybrids Pelophylaxesculentus. P. lessonae and P. esculentus 
form stable complexes (L-E complexes) in which P. esculentus are hemiclonal. In L-E complexes all the transmitted 
genomes by P. esculentus carry deleterious mutations which are lethal in homozygosity. 

Results: We analyze, by means of an individual based computational model, L-E complexes. The results of simulations 
based on the model show that, by eliminating deleterious mutations, L-E complexes collapse. In addition, simulations 
show that particular female preferences can contribute to the diffusion of deleterious mutations among all P. 
esculentus frogs. Finally, simulations show how L-E complexes react to the introduction of translocated P. ridibundus. 

Conclusions: The conclusions are the following: (i) deleterious mutations (combined with sexual preferences) 
strongly contribute to the stability of L-E complexes; (ii) female sexual choice can contribute to the diffusion of 
deleterious mutations; and (iii) the introduction of P. ridibundus can destabilize L-E complexes. 
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Background 

Lake frog (Pelophylax ridibundus Pallas, 1771) and pool 
frog [Pelophylax lessonae Camerano, 1882) can mate 
producing the hybrid edible frog {Pelophylax esculentus 
Linneus, 1758). P. esculentus can coexist with one or both 
of the parental species giving rise to mixed populations. 
Usually, the genotypes of P. ridibundus, P. lessonae and P. 
esculentus are indicated by RR, LL, and LR, respectively. 
In Europe there are mixed populations containing P. ridi- 
bundus and P esculentus individuals, called R-E systems, 
populations with P lessonae and P esculentus individu- 
als, called L-E systems, and populations with all three 
species. Due to the eastern origin of P ridibundus, R-E 
complexes are frequently found in Eastern Europe, while 
L-E systems are widespread throughout the rest of Europe 



"Correspondence: barbuti@di.unipi.it 

1 Dipartimento di Informatica, Universita di Pisa, Largo B. Pontecorvo, 3, 561 27 
Pisa, Italy 

2 Museo di Storia Naturale, Universita di Pisa, Via Roma, 79, 5601 1 Calci (Pisa), 
Italy 



[1-5]. Hybrids in these populations reproduce in a par- 
ticular way, called hybridogenesis [1,6-12]. Hybridogenesis 
consists of a gametogenetic process in which the hybrids 
exclude one of their parental genomes premeiotically, and 
transmit the other genome, clonally, to eggs and sperm. 
For example, in L-E complexes, P esculentus hybrids have 
both the genomes of the parental species, L and R, but they 
produce only R gametes. 

This mode of reproduction requires hybrids to live sym- 
patrically with the parental species whose genome has 
been eliminated. In this way hybrids in a L-E system elim- 
inate the L genome thus producing P esculentus when 
mating with P lessonae, and generating P ridibundus 
when mating with other hybrids. Usually P ridibundus 
generated in L-E complexes are inviable due to deleteri- 
ous mutations accumulating in the clonally transmitted R 
genome [13-17]. Analogously, in R-E systems there is a 
tendency during hybrid gametogenesis to eliminate the R 
genome; as with L-E systems, P lessonae, the offspring of 
hybrids, are often inviable. 
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In natural L-E complexes, the inviability of offspring 
of P. esculentus x P. esculentus matings is evidenced 
by the absence of adults P. ridibundus. Experimental 
crosses between coexisting hybrids (from localities sam- 
pled throughout the range of L-E populations) also show 
that such offspring are inviable [6,14,16,17]. These stud- 
ies have also revealed that the same hybrid individuals, 
producing inviable progeny, produce viable progeny when 
crossed with either parental species or with hybrids from 
different regions. The lethality of natural hybrid x hybrid 
matings is thus neither the result of hybrid sterility nor 
the inherent consequence of the hemiclonal reproduc- 
tive mode. Guex et al. present two simple hypotheses, 
both explaining the observed inviability of P. esculentus x 
P esculentus progeny by the load of deleterious muta- 
tions on the clonally transmitted R genomes [16]. Quoting 
from their paper, the hypotheses are "(1) inviability is 
caused by homozygosity for recessive deleterious muta- 
tions at particular gene loci; or (2) inviability is caused 
by a general deterioration of non-recombining R genomes 
through Mullers ratchet, reflecting different hemiclone- 
specific sets of incompletely recessive mutations, which 
leads to lethality when two such deteriorated R genomes 
are combined". Their conclusion is that the hypotheses 
are not mutually exclusive, however there is evidence to 
support the plausibility of the first hypothesis: Mullers 
ratchet generates deleterious mutations in relatively ran- 
dom places in the genome, which are then likely to be 
different in different geographical areas. The study in 
[16] suggests that, in some cases, single lethal mutations 
may have the same effect as the accumulation of deleteri- 
ous ones. However, most studies on P esculentus fitness, 
suggest that when in a heterozygous state the effect of 
deleterious mutation is not significant. 

Due to the inviability of P esculentus x P esculen- 
tus offspring, P esculentus populations cannot survive 
alone, but must act as a sexual parasite of one of the 
parental species. This dependency can be avoided only by 
all-hybrid populations in which the presence of triploid 
and tetraploid individuals leads to recombination among 
homolog parental chromosomes [18-23]. This recombina- 
tion is able to purge, at least partially, deleterious muta- 
tions from genomes, thus producing viable offspring of 
hybrids. Due to its wide distribution, we will consider L-E 
complexes, in which P esculentus and P lessonae coexist. 
In such a complex, the reproductive pattern is shown in 
Table 1, where the subscript y indicates the male sexual 
chromosome. 



Table 1 Reproductive pattern of water frogs 





LL 


LR 


LyL 


LyL LL 


L y R LR 


LyR 


LR 


RR not viable 



The Y chromosome determines the sex of frog males 
and can occur only in the L genome, due to primary 
hybridization which involves, due to size constraints, P 
lessonae males and P ridibundus females. Only one of 
the three possible matings resulting in viable offspring 
produces LL genotypes (Table 1). This would give an 
advantage to P esculentus which could outnumber P 
lessonae and eventually eliminate them. This situation 
would also eventually result in an extinction of P escu- 
lentus which cannot survive without the parental species. 
In addition to their relative abundance which is pro- 
moted by the above reproductive pattern, P esculentus 
have other advantages. Although in many cases they show 
either no differences or intermediate characteristics com- 
pared to their parental species [24-27], P esculentus show 
behavioural differences [28,29], and, have, by heterosis, a 
greater fitness than the parental species in certain aspects 
[13,15,30-36]. The combination of the relative abundance 
and heterosis should out-compete P lessonae in L-E com- 
plexes. The widespread distribution of L-E complexes, 
although with different percentages of hybrids, reveals 
that there are mechanisms which contribute to the sta- 
bility of such complexes [37-40]. Of these mechanisms, 
sexual selection seems to be one of the most important. 
In fact, P esculentus females prefer (either overtly or cryp- 
tically) P lessonae males than males of their own species 
[41-45]. Many mathematical and computational models 
have studied the influence of sexual selection in the evolu- 
tion of populations [46-54]. In addition some models have 
focused on sexual selection in complexes in which some 
form of clonal reproduction exists [55-57]. The models 
in [57-59] show how female preference is able to stabilize 
L-E complexes by counterbalancing both heterosis and 
the reproductive advantage of £ esculentus. Other factors, 
such as reproductive performance, in conjunction with 
sexual choice can increase the stability of L-E complexes 
[60]. 

Using an individual-based computational model, in this 
paper we study three problems. The first is how deleteri- 
ous mutations contribute to the stability of L-E complexes. 
The second concerns how, in L-E complexes, deleterious 
mutations can diffuse in the R genomes of the whole P 
esculentus population. The third is the invasiveness of P 
ridibundus in L-E complexes. 

Regarding the first problem, the aim is to investigate 
whether deleterious mutations on the R genome con- 
tribute, together with female preferences, to the stability 
of L-E complexes. 

As for the diffusion of deleterious mutations in the pop- 
ulation, an interesting hypothesis is proposed in [61], in 
which deleterious mutations can influence female pref- 
erences. From the literature we know that females of P 
esculentus have a strong preference for P lessonae males. 
Vorburger et al. in [61] suggest that, of £ esculentus males, 
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P esculentus females may prefer those with mutations on 
the R genome. Such mutations could make the affected 
loci on the R genome dysfunctional, thus producing a 
more u lessonae Aike" genotype. 

Finally, regarding the third problem, Vorburger and 
Reyer in [9] suggest that the introduction of P ridibun- 
dus can either provoke the collapse of L-E populations, or 
result in a replacement by P ridibundus of both P lessonae 
and P esculentus, leading to a mono-specific population. 

In order to gain further insight into the three problems, 
we attempt to answer three questions. 

- Is the role of deleterious mutations necessary for the 
stability of L-E complexes? 

- How can a stable L-E complex be obtained? 

- What is the effect of introducing P. ridibundus into 
L-E complexes? 

Methods 

The model 

To study the interaction between populations of P 
lessonae, P esculentus and P ridibundus we developed an 
individual-based model. To answer the three questions 
above, we started with a simple model (for the first ques- 
tion) and then extended it step by step (to tackle the 
second and third questions). 

In the simplest model, we consider diploid individuals 
whose genotype is represented by two chromosome types: 
L and R. Chromosomes R can contain deleterious muta- 
tions (represented by R d ), and only chromosomes L can 
have the sex-determining chromosome Y (represented by 
L y ). Thus the possible genotypes are: LL, L y L, LR, L y R, LR d , 
L y R d ,RR,RR d ,R d R d . 

The fitness of a genotype g, T{g), is computed as fol- 
lows: 



0 if g = R d R d 

_( a-c(g)) 2 \ 



otherwise 



where a is a parameter measuring the strength of the 
ecological selection (smaller values of a correspond to a 
stronger selection). In the simulations we use two differ- 
ent values for a, a = 0.4, which corresponds to a hard 
environment, and o = 0.6, which corresponds to a weaker 
selection. 

Function c(g), is the chromosomes fitness, defined as 
follows: 



c(g) = 



1 if ge{LR, L y R, LR d ,L y R d ] 

l — Sh tig e {LL, LyL} 
l-8 h -8 e ifge{RR,RR d } 



where 8^ and 8 e describe the fitness decrement associated 
with homozygous genotypes (which do not gain from het- 
erosis), and P ridibundus genotypes, respectively. The use 
of a further decrement, 8 e , in the chromosome fitness of 



P ridibundus, derives from the fact that L-E complexes 
usually live in pools and marshes, where P ridibundus 
are less fit. In the following sections, we use 8^ = 0.2 
and 8 e will assume the values 0.0, 0.2, and 0.4. 8 e = 
0.0 means that P ridibundus have the same fitness as P 
lessonae (i.e. the environment includes niches for both 
species), while 8 e = 0.4 represents the fact that P ridibun- 
dus are strongly disadvantaged compared to P lessonae 
(i.e. the environment consists of a typical P lessonae 
habitat). 

We consider that the populations have a reproductive 
season each year. During this season all the females mate. 
Female sexual choice is implemented by a best-of-n selec- 
tion procedure of males [52,62], i.e. a female mates with 
the most preferred of n randomly chosen males in the 
population. The best-of-n procedure, which is a usual 
computational method to take into account female prefer- 
ences, can be described as follows. Consider a population 
of males of k different types, m\, m^, . . . m^, with num- 
bers of individuals N mi ,N m2 , . . .N mk , respectively, and a 



N mi + N m7 + . . . + N, 



V W2 



m k - 



total male population N m 
Consider a female of type / with different preference val- 

f f f 
ues for k kinds of males, p mi>P m 2 > — - >P m k > such that 

f f f 

Pm\ < Pm 2 < . . . < Pm k < According to the best-of-n 

procedure, the female chooses the best (following her 
preferences) of the n ones randomly taken from the pop- 
ulation of all males. Thus the probability, prob m , of a 
female of type/, choosing a male of type mu i e [0,k], is 
given by the probability of randomly taking no mj males, 
f f 

with p m . > p m , times the probability of having at least 
one mt male of the taken ones. Using the hypergeometric 
probability distribution, pro\J m is given by the following 
formula. 



prov mi 



bL; = 



(N m ,+N rm +...+N rn 



i N n) 



*m 2 - 

n—i 



■N mi +N m2 +...+N m 



) 



Note that the greater n is, the greater the strength of 
the female preference. Increasing the n value will lead 
to a female choosing from a greater number of males, 
thus mimicking the behaviour of a more discriminating 
female. In order to obtain stable complexes, we assume 
a species-specific female preference, in particular, follow- 
ing the studies in [59,60], we assign a stronger prefer- 
ence to females of the parental species than the hybrid 
females. P esculentus females have the same behaviour as 
P lessonae, thus competing for the same kind of males. 
Likewise, P ridibundus females prefer more "ridibundus- 
like" males. Hereafter we call u lessonae preference" this 
kind of female preference, because P lessonae males are 
the most preferred. The number of n candidates, which 
each female chooses from, is set to 30, 15, and 30 for 
the P lessonae, P esculentus, and P ridibundus females, 
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respectively. Choosing from among 30 candidates has a 
biological flavour. If males are distributed in the envi- 
ronment with a density of one male per square meter, 
each female must swim in a circle of 10 meters of diam- 
eter in order to check out 30 males, a distance which 
is reasonable for a frog. In [57-59] it is shown that sta- 
ble populations are only found when the preference of P 
lessonae females is greater than the preference of P escu- 
lentus females. In these papers it is shown that, under the 
above assumption, many different values of female prefer- 
ence lead to stable complexes. We have performed many 
simulations by leaving the value of n for P. lessonae and P. 
ridibundus unmodified, and by varying the value of n for 
P esculentus. We have found that if the ratio of the value of 
P esculentus to the value of the parental species belongs to 
the interval [0.03, 0.7] we obtain stable complexes, when 
P ridibundus are inviable. This result is analogous to the 
ones in [57-59]. Because different values in the interval 
[0.03, 0.7] affect only the percentage of hybrid frogs in the 
final stable population, the choice of n does not change the 
overall dynamics of the population. For this reason we use 
the non-extreme value 0.5. 

Offsprings genotypes are obtained from the gamete 
combination of the parents. 

The reproductive season is followed by a viability selec- 
tion. During this phase the probability that an individual 
of genotype g survives, p SU rv(g)> is given by a slight modi- 
fication of the Beverton-Holt model: [63-68]: 



where b is the average number of offspring for females 
that can reach the stage of adults, 0 is the percentage of 
females in the population, N is the number of individ- 
uals in the population competing for the resources, and 
JC(g) is the carrying capacity associated with the genotype 
g. JC(g) is given by F(g)K,$ where /Co is the maximum 
carrying capacity of the environment. In all our simula- 
tions we assume b = 6 and /Co = 3000. The standard 
Beverton-Holt model is modified due to the fact that we 
consider overlapping generations and we apply the via- 
bility selection based on survival probability, not only to 
young tadpoles but to all the individuals in the population. 
The simplest model is used to answer the first question 
by performing both simulations in which all R genomes 
carry deleterious mutations and simulations in which all R 
genomes are free from deleterious mutations. 

To answer the second question, "How can a stable L-E 
complex be obtained?" we simulate the diffusion of dele- 
terious mutations in the population, starting with an L-E 
complex, composed by P lessonae, in which there are 
only a few P esculentus individuals (without mutations on 
the R genome). We consider a mutation rate, /x, which 



gives the probability of adding a new deleterious muta- 
tion on the R genome of an offspring (in the simulations 
[A is set both to 10 -4 and to 10 -5 ). We consider differ- 
ent "stages" in the accumulation of mutations in the R 
genome. For the sake of simplicity we consider only three 
stages: R, R d i, and R d . R is the genome without muta- 
tions, Rd\ is the genome with a non lethal accumulation 
of mutations, and, as before, R d is the final stage of accu- 
mulation (lethal in homozygous individuals). With this 
scenario we have P ridibundus females with the follow- 
ing possible genotypes (in order of decreasing fitness): 
RR, RR dl , RR d , R d iRdi> Rdi*d> Rd*d> where R d R d is the 
only lethal genotype. The fitness decrease of each geno- 
type with respect to the previous one is given by 8 m . In 
the simulations we set 8 m = 0.04. In this scenario, fol- 
lowing [61], we assume that P lessonae and P ridibundus 
females have a strong preference for males of their own 
species, while of P esculentus males, females prefer males 
with a more "lessonae-like" genotype (males with muta- 
tions on the R genotype). We call this kind of female 
preference "lessonae-like preference". We have P esculen- 
tus males with the following genotype: L y R, L y R d i, L y R d , 
which, according to the "lessonae-like preference" are in 
order of increasing preference by the female. 

Finally, to answer the third question, "What is the effect 
of introducing P ridibundus into L-E complexes?" we need 
to simulate the effect of the introduction of P ridibundus 
males and females in an L-E complex. In these simulations 
we consider both R genomes without mutations; this is 
because we assume the absence of the Muller ratchet in 
the (sexually reproducing) introduced P ridibundus and 
R d genomes generated by the stable L-E complex. The 
further extension we consider is the possibility of having 
P ridibundus males, i.e. the possibility of having the Y 
chromosome on R genomes as well: R y and R yd . 

We introduced a limit for the lifespan of individuals, 
thus all individuals exceeding 10 years of age are removed 
from the population. Removing old frogs from the system 
avoids that extremely fit individuals survive indefinitely 
because the viability selection is not able to remove them. 
This situation can easily happen in models with overlap- 
ping generations which do not consider deaths due to 
aging. 

The parameters used in the model, with their meaning 
and values, are reported in Table 2. 

Implementation 

The process was formalize with CLS [69-75] and then 
implemented in C++ language using Microsoft Visual 
Studio 2010. Each simulation runs for 2000 iterations 
(more when necessary) with a carrying capacity of 3000 
individuals. A single simulation takes nearly lOsec on a 
worstation equipped with Intel i5 3.2Ghz processor and 
8GB RAM. 
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Table 2 Parameters used in the model 



o Selection strength 0.4 - 0.6 

Sh Fitness decrease due to homozygosity 0.2 

S e Fitness decrease for environmental reasons 0.0 - 0.4 

8 m Fitness decrease for each step of mutation accumulation 0.04 

n Number of attempts per female in choosing males, P. lessonae, P. esculentus, P. ridibundus 30, 15,30 

b Number of surviving tadpoles for each female 6 

/Co Environmental carrying capacity 3000 

li mutation rate 10 -5 , 10 -4 



Results 

Is the role of deleterious mutations crucial for the stability 
of L-E complexes? 

The first step in studying the role of deleterious muta- 
tions for the stability of L-E complexes is to consider an 
initial population composed of P. lessonae frogs and small 
percentages of P. esculentus: 5%, 10%, and 20%, respec- 
tively. Note that an initial situation with a large number 
of P. lessonae individuals favours the stabilization of the 
complex (there is no possibility of an early collapse due 
to the greater fitness of P esculentus). We assume that 
all P esculentus individuals carry the deleterious muta- 
tions on the R genome, that is P ridibundus females are 
not viable and they do not appear in the population. We 
assume the "lessonae preference" for females, based on 
the values 30, 15, and 30 for the best-of-n procedure for 
the three species, as described in the previous section. 
Finally we consider 8^ = 0.2. We perform simulations with 
a = 0.4 and a = 0.6. For each combination of parame- 
ters o and 8 e we performed 100 simulations, the possible 
outcomes of which are either a stable L-E complex or 
the collapse of the whole population. In all the simula- 
tions the populations evolve towards a stable L-E com- 
plex, following a typical population composition pattern 
(Figure 1). 

The result of these simulations is not surprising. Essen- 
tially the results are in accordance with those in [58,59], 
showing that female preference is a strong stabilizing force 
for L-E complexes. 



In order to investigate the role of deleterious mutations 
in the stability of L-E complexes, we consider the same 
parameters as the previous simulations, but we remove 
the deleterious mutations from all the R genomes of the 
initial population. In addition, we set 8 e to 0.2 and 0.4, that 
is we consider that all P ridibundus born from hybrids 
are always disadvantaged compared to P lessonae and P 
esculentus. We also set the mutation rate /x equal to 0, in 
order to prevent deleterious mutations, considering only 
a "mutation-free" population. We can observe that with 
all the parameter combinations and the initial percent- 
ages of P. esculentus, in all the simulations the population 
eventually collapses. If viable P. ridibundus females are 
produced, the reproductive pattern becomes as depicted 
in Table 3. The table highlights that this reproductive pat- 
tern generates a numerical disadvantage for P. lessonae, 
the population of which decreases. The decrease in the P. 
lessonae population has, as a consequence, a decrease in 
produced L gametes, which, in turn, results in a higher 
production of P. ridibundus. Thus the population of P. 
ridibundus females grows and eventually out-competes 
the other species, despite the weaker fitness of P. ridibun- 
dus females compared both to P. esculentus (8^ = 0.2) and 
to P. lessonae (8 e = 0.2 and 8 e = 0.4), see Figure 2. Of 
course a population of only P. ridibundus females cannot 
survive. 

Note that the higher the percentage of P. esculentus in 
the initial population, the smaller the number of genera- 
tions before the collapse. By varying either the strength of 
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Table 3 Reproductive pattern of water frogs without 
deleterious mutations 





LL 


LR 


RR 


LyL 


L y L LL 


L y R LR 


L y R LR 


LyR 


LR 


RR 


RR 



female preferences or the strength of ecological selection, 
we obtain different values for the number of generations 
before the collapse, but an evolution towards collapse 
remains the overall trend of the system. Even when P. 
ridibundus are at a serious disadvantage, 8 e = 0.4, RR 
females are able to compromise the stability of L-E com- 
plexes. Of course, we amplify this effect by decreasing 
8 e , thus we do not show the population dynamics with 
8 e = 0.0. In this case, the point of collapse is reached very 
quickly. 

How can a stable L-E complex be obtained? 

In this section we study the effect of both "lessonae pref- 
erence" and "lessonae-like preference" of R esculentus and 
R lessonae females in the diffusion of deleterious muta- 
tions. We set the mutation rate, /x, either to 10 -4 or 
to 10 -5 . We consider three possible stages of deleteri- 
ous mutation accumulation in each R genome (R, Rji, 
and Rd), thus any mutation event determines the pas- 
sage from one stage to the next. We start the simulations 
with two different initial populations. The consistency 
of the initial R lessonae population is the same in all 
the simulations, 2700 individuals, but the number of R 
esculentus individuals is set, initially, either to 10 or 100 
individuals. We also perform simulations with two values 
for 8 et which lead to a decrease in fitness of mutation- 
free R ridibundus compared to R lessonae. We set 8 e to 
either 0.2 or 0.4, i.e. we consider that the environment is 
either weakly penalizing for R ridibundus or is a typical 
R lessonae habitat, which is not suitable for R ridibundus 
frogs. In addition, any further mutation accumulation on 
the R genomes of R ridibundus decreases their fitness by 
8 m = 0.4. 



Firstly, we consider the "lessonae preference". Accord- 
ing to this preference pattern, R esculentus and R lessonae 
females prefer R lessonae males, and, do not choose 
among R esculentus males. The results of the simulations 
show that in our model both with 10 and with 100 R escu- 
lentus mutation-free individuals in the initial population, 
and with any value for 8 e and cr, the lethal deleterious 
mutation diffuses slowly. This slow diffusion is not able to 
prevent the production of a sufficient number of viable R 
ridibundus females, which leads to the collapse of the pop- 
ulation with all possible parameter values. Figure 3 shows 
a typical dynamic of the population while Figure 4 shows 
the diffusion of mutations in the R genome of R esculentus. 
The results of the simulations with a = 0.4 and fi = 10 -4 
are shown in Table 4. 

Secondly, we consider the "lessonae-like preference". 
According to this preference pattern R esculentus and 
R lessonae females prefer R lessonae males, and, of R 
esculentus males, they prefer those with a greater accumu- 
lation of deleterious mutations. The results of simulations 
confirm that female choice can lead to the diffusion of 
deleterious mutations in the R esculentus genomes if those 
mutations silence R ridibundus traits and result in more 
"lessonae-like" R esculentus males. In this scenario, sta- 
ble complexes can be obtained with 8 e = 0.4, while a 
higher fitness of R ridibundus leads to the extinction of 
the population (Table 5). Starting with mutation-free R 
esculentus individuals, the diffusion of mutations passes 
through different stages of mutation accumulations. The 
pattern of mutation diffusion is that mutation rates con- 
tribute to deleterious mutations, while the "lessonae-like 
preference" favours the production of offspring of R escu- 
lentus males with a stronger mutation accumulation on 
their R genome. In the progression of mutation accumu- 
lation, RR genotypes gets more and more inviable while 
LR males become more and more attractive. When 8 e = 
0.4, Figure 5 shows a typical dynamic of the popula- 
tion, while Figure 6 shows the diffusion of mutations in 
the R genome of R esculentus. The results of the sim- 
ulations with a = 0.4 and fi = 10 -4 are shown in 
Table 5. 



3000 
2500 




Figure 2 Results of typical simulations with an initial population composed of P. lessonae frogs and a percentage of 5% mutation-free P. 
esculentus. a = 0.4, 8 e = 0.2, \x = 0. 
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What is the effect of the introduction of P. ridibundus in L-E 
complexes? 

In this section we analyze the effect of introducing P ridi- 
bundus into an L-E complex. This scenario really happens 
in natural environments due to the importation in Western 
Europe of P ridibundus for commercial purposes. 

In order to study the effect of the translocation of P. ridi- 
bundus, we performed simulations by varying the fitness 
of the introduced frogs (8 e = 0.0, 0.2, 0.4). 8 e = 0.0 means 
that the environment does not put P ridibundus at a dis- 
advantage with respect to P lessonae. In these simulations 
we consider a strong selection strength, a = 0.4. 

We study the effect of introducing, in a stable L-E com- 
plex, a percentage of both 5% and 10% of P ridibundus, 
males and females. We consider three different situations: 
i) P ridibundus are introduced into the typical environ- 
ment for P lessonae (8 e = 0.4), ii) P ridibundus are intro- 
duced into a mixed (lake/marsh) environment (8 e = 0.2), 
and iii) P ridibundus are released into an environment that 
is fit for them (possibly close to resident frogs) (8 e = 0.0). 
We have four possible outcomes: stable L-E systems, sta- 
ble L-E-R systems, stable P ridibundus populations, and 
the collapse of the population (Table 6, Figures 7, 8, and 
9). In the figures we show the dynamics for a = 0.4. For 
o = 0.6 we obtain analogous evolutions. 

If the fitness of P ridibundus is equal to the fitness of P 
lessonae (8 e = 0.0), in most cases the population becomes 



a mono-specific population of P ridibundus. Note that 
there are many cases in which the three species coexist 
at the end of the simulations, which we discuss explic- 
itly in the following section. If the fitness of P ridibundus 
decreases (8 e = 0.2), the introduced frogs do not survive 
for long. However, before their extinction P ridibundus 
can mate with P lessonae and P esculentus, thus introduc- 
ing mutation-free R genomes into the hybrid population. 
At this point, viable P ridibundus females born from 
matings between hybrids lead the population to collapse. 
Finally, when the fitness of P ridibundus is very low (8 e = 
0.4), P ridibundus frogs are immediately expelled from the 
system. 

Discussion 

From the results of the previous section we can deduce 
the overall dynamics of L-E complexes with viable P 
ridibundus females. In general, the presence of viable P 
ridibundus females significantly changes the reproductive 
outcome of L-E complexes, and the generation of off- 
spring becomes as depicted in Table 3. For the sake of 
simplicity let us consider hypothetical complexes in which 
the number of females is the same for the three species. 
For such populations the production of P lessonae off- 
spring passes from 33.33% in stable L-E complexes (recall 
that in such complexes P ridibundus are inviable and do 
not survive) to 16.66% in complexes where P ridibundus 
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Figure 4 Results of typical simulations with an initial population composed of P. lessonae frogs and 10 mutation-free P. esculentus. In the 

simulation we assume a "lessonae preference" for P. esculentus and P. lessonae females, o = 0.4, 8 e = 0.4, fi = 1 0 -4 . R, R^d, and Rd are the number of 
R genomes in the P. esculentus population with three different levels of mutation accumulations. 
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Table 4 Outcomes of 100 simulations with an initial P. 
esculentus population of 1 0 and 1 00 individuals without 
deleterious mutations on the R genome 



a) 8 e = 0.2 


Initial P. esculentus pop. 




10 


100 


Stable L-E system 




0 


0 


LL population 




0 


0 


Collapse 




100 


100 


b) S e = 0.4 


Initial P. esculentus pop. 




10 


100 


Stable L-E system 




0 


0 


LL population 




0 


0 


Collapse 




100 


100 



In the simulations we assume a "lessonae preference" for P. esculentus and P. 
lessonae females. The parameters are: a = 0.4, \i = 1 0" 4 . 



survive. Thus the viability of P. ridibundus decreases the 
relative abundance of P. lessonae offspring. The decrease 
in the relative production of P lessonae offspring leads to 
a decrease in P lessonae adults in the future. This, in turn, 
causes a decrease in the production of L gametes, which 
are produced only by P lessonae frogs. This process, if 
not stopped by some external reasons, results in a trend 
towards the extinction of P lessonae, A population com- 
posed only of P esculentus individuals and P ridibundus 
females cannot survive, because no gametes with the Y 
chromosome can be generated. 

The discussion which follows, regarding the three ques- 
tions mentioned previously, is based on the trend towards 
extinction mentioned above, modulated for example by 
the accumulation of mutations, or the introduction of 
viable P ridibundus males. 



Table 5 Outcomes of 100 simulations with an initial P. 
esculentus population of 1 0 and 1 00 individuals without 
deleterious mutations on the R genome 



a) 8 e = 0.2 


Initial P. esculentus pop. 




10 


100 


Stable L-E system 




0 


0 


LL population 




0 


0 


Collapse 




100 


100 


b) S e = 0.4 


Initial P. esculentus pop. 




10 


100 


Stable L-E system 




76 


73 


LL population 




1 


0 


Collapse 




23 


27 



In the simulations we assume a "lessonae-Wke preference" for P. esculentus and P. 
lessonae females. The parameters are: o = 0.4, /x = 1 0" 4 . 



Deleterious mutations are necessary for the stability of L-E 
complexes 

We showed in Section Is the role of deleterious mutations 
crucial for the stability of L-E complexes?' that by using 
the "lessonae preference" we essentially obtain the same 
results as [58,59]. If there are no viable P ridibundus off- 
spring, the system evolves towards stability regardless of 
the strength of the selection. The same complexes, with 
the same values of female preferences, but without delete- 
rious mutations on R genomes will collapse irrespectively 
both of the selection strength and the fitness of viable P 
ridibundus. In the dynamics of the population towards 
collapse, there are roughly three phases. In the first phase, 
the effect of sexual selection and the abundance of P 
lessonae males maintains a low number of P ridibundus 
offspring. Although no mutations are present in the P 
esculentus genome, a very few P ridibundus females are 
generated because mating between hybrids is rare, due to 
female preferences. In the second phase, the number of 
hybrids increases because of their fitness as a result of het- 
erosis, while the number of P lessonae decreases. In this 
phase the increased number of hybrids facilitates mating 
which thereby produces viable P ridibundus females. In 
the third phase, P ridibundus females act as sexual par- 
asites for hybrids, and their number grows because of 
P esculentus x P esculentus and P esculentus x P ridi- 
bundus matings. Note that these matings produce only 
P ridibundus females. Viable P ridibundus females will 
mate, preferably, with P esculentus males (which have a 
more "ridibundus phenotype") producing a bigger popu- 
lation of ridibundus females for any further generation. In 
the model we do not consider male frogs as a limiting fac- 
tor for reproduction, that is one male is able to fertilize the 
eggs of a unlimited number of females. We consider also 
that males will not make any specific choice of females 
[45]. The number of P esculentus decreases because of 
the reduced number of P lessonae. This phase ends with 
a population of all P ridibundus females which quickly 
collapses (Figure 2). 

Consequently, if stability is maintained by female pref- 
erences, collapse only can be prevented if the initial popu- 
lation of P esculentus is affected by deleterious mutations 
on the R genomes, which prevents the birth of viable P 
ridibundus. Thus we can conclude that, in a system in 
which the parameters used for the stability are essentially 
female preferences and viability of P ridibundus, dele- 
terious mutations are necessary for the stability of the 
complex. These results help to clarify why natural L-E 
complexes generating viable P. ridibundus are extremely 
rare, and in most cases the percentage of viable P. ridibun- 
dus is not significant [15]. It is difficult (perhaps impos- 
sible) for L-E complexes to persist if they generate viable 
P ridibundus. Starting with an initial population in which 
all P. esculentus frogs carry deleterious mutations, the 
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complex evolves towards a stable configuration (Figure 1). 
The possible collapse of the population occurs only if P 
esculentus out-compete P lessonae, however this evolu- 
tion seldom takes place because of the large number of P. 
lessonae in the initial population. The results show that, 
with the same selection strength, L-E complexes evolve 
towards the same percentages of the two frog species, 
whatever the initial percentages. 

Our results regarding the stability of L-E complexes 
both with deleterious mutations in all R genomes and with 
a "lessonae preference" are similar to those in [57-59]. We 
show that deleterious mutations strongly influence the 
stability of L-E complexes. Thus deleterious mutations are 
not only a secondary consequence of Mullers ratchet but 
they have an important role in the stability of complexes. 
In this paper we highlight that, if we consider both hetero- 
sis and females preferences as the only stabilizing forces 
of L-E complexes, the lack of deleterious mutations drives 
such populations to collapse. Thus neither female prefer- 
ences nor deleterious mutations are sufficient to maintain 
the stability of L-E complexes however, in this scenario, 
each is necessary for stability. 

Female preference can contribute to obtaining stable L-E 
complexes 

The above discussion highlights the important role of 
deleterious mutations in stable L-E complexes. But how 
can stable L-E complexes be obtained? 



Our results show that, in order to reach stable L-E 
populations, there must be forces that drive the diffu- 
sion of deleterious mutations. Even with a fast mutation 
rate, 10 -4 , Mullers ratchet alone is not sufficient for dif- 
fusing the deleterious mutations in the population. The 
"lessonae preference" does not allow P lessonae and P 
esculentus females to choose from among P esculentus 
males. Thus the diffusion of deleterious mutation among 
P esculentus individuals and P ridibundus females is not 
guided by female preferences, but only by the mutation 
rate (essentially by Mullers ratchet). In all our simula- 
tions this diffusion turns out to be very slow, thus viable P 
ridibundus females are generated before Mullers ratchet 
accumulates lethal mutations on all R genomes. This 
intermediate phase, with a sufficient number of viable P 
ridibundus females, is responsible for the collapse of the 
whole system. 

Selecting R genomes with a higher mutation accu- 
mulation is accelerated by the "lessonae-\ike preference". 
In this case the production of offspring with greater 
mutation accumulation on the R genome is favoured, 
and consequently the production of fit P ridibundus 
females is lowered. To understand this process we need 
to consider that mutation accumulation on R genomes 
decreases the fitness of P ridibundus females, but it does 
not affect the fitness of P esculentus in which dysfunc- 
tional R genomes are counterbalanced by "healthy" L 
genomes. 




Figure 6 Results of typical simulations with an initial population composed of P. lessonae frogs and 10 mutation-free P. esculentus. In the 

simulation we assume a "lessonae- 1 ike preference" for P. esculentus and P. lessonae females, o = 0.4, 8 e = 0.4, fi = 1 0 -4 . R, R^d, and Rj are the 
number of R genomes in the P. esculentus population with three different levels of mutation accumulations. 
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Table 6 Outcomes of 100 simulations in which in a stable 
L-E complex a percentage of 5% of P. Hdibundus, males 
and females, are introduced 





8 e = 0.0 


8 e = 0.2 


8 e = 0.4 


Stable L-E system 


0 


46 


100 


Stable L-E-R system 


32 


0 


0 


Stable P. ridibundus population 


39 


0 


0 


Collapse 


29 


54 


0 



The three columns correspond to different values of 8 e , 0.0, 0.2, and 0.4, 
respectively, a = 0.4, \x = 0. 



Another important point is that a significant parameter 
in the diffusion of deleterious mutations is the fitness of P 
ridibundus to the environment. If this fitness is too high, 
too many viable P ridibundus females are produced before 
significant mutation accumulations. We know that a high 
number of such viable females will lead the population to 
collapse. 

For computational purposes we have considered only 
three stages of mutation accumulation on the R genomes. 
This is an approximation of the Mullers ratchet which, in 
most cases, operates through a huge number of mutations. 
We approximate the Mullers ratchet by decreasing the 
mutation rate so that a mutation in the model corresponds 
to many mutations in real genomes. Following the estima- 
tion in [76], we assume that, in an eukaryote organism, the 
mutation rate in a whole genome during sexual reproduc- 
tion is in the interval [ 3 x 10 -2 , 9 x 10 -1 ]. Many of these 
mutations are either not significant or not deleterious. 
Values of the mutation rate in the interval [ 10 -5 , 10 -4 ], 
used in our model, take both the above considerations into 
account. 

Our study differs significantly with the results of other 
authors with regard to the diffusion of mutations. The 
models in [58-60] provide an extensive insight into the 
reasons for the stability of L-E complexes starting from 
the assumption that deleterious mutations are present in 
all the R genomes in the population. Our model builds 
on the previous ones by assuming sexual choices in the 



populations. However, it differs by considering a popula- 
tion in which deleterious mutations are not present, but 
are generated according to a mutation rate. In addition, 
only when this accumulation reaches a given threshold 
does it become lethal. This leads us to conclude that sexual 
selection not only stabilizes the complexes, but can also 
force mutation diffusion. 

Note that our simulations do not enable us to prove 
the hypothesis suggested in [61]. Computational and 
mathematical models, without subsequent experimen- 
tal support, can only be used to rule out incorrect 
hypotheses - they cannot prove correct ones. Compu- 
tational and mathematical models can only state that a 
hypothesis is plausible. In the case of L-E complexes, the 
stabilization period is so long that no real experiment 
can support a hypothesis on its stabilization, however 
our model suggests that the hypothesis in [61] regard- 
ing female preference, could plausibly lead to stable L-E 
complexes. 

Invasion of translocated P. ridibundus 

Another point that we study with our model is the con- 
sequence of introducing P. ridibundus into stable L-E 
complexes. P. ridibundus can mate both with P. esculen- 
tus, producing P. ridibundus, and with P. lessonae (primary 
hybridization), producing P. esculentus. Primary hybrids 
can have low fertility rates [77], thus their contribution 
to the dynamics of the population is low. In our model 
we take account of this low contribution by decreasing 
the possibility of producing primary hybrids. This is done 
by setting the female preferences of P. ridibundus in a 
way that P lessonae males are seldom chosen (the value 
of n, in the best-of-n procedure, is set to 30). On the 
other hand, the preference of P lessonae females is mainly 
for males of their own species. We assign the same fit- 
ness both to the introduced P ridibundus and to those 
generated by matings of P ridibundus with P esculentus. 
This is in accordance with the semi-natural experiments 
in [17]. However, in some simulations we use a P ridi- 
bundus fitness, which is lower than the fitness of the 
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resident P lessonae and P esculentus because we consider 
that a marshy environment with a low oxygen level, where 
P. lessonae and P esculentus live, is less suitable for P 
ridibundus. 

The results show that, as predicted in [9], the intro- 
duced P ridibundus often out-compete the other species 
resulting in a mono-specific population, when their fitness 
is comparable to the resident populations. Although the 
introduction of P ridibundus results in new R hemiclones, 
which contribute to the genetic diversity of hybrids, our 
results do not support the hypothesis presented in [5] that 
this genetic diversity can stabilize the hybridogenetic sys- 
tem. If the fitness of P ridibundus is competitive with the 
fitness of the resident population (8 e = 0.0), P ridibun- 
dus males will survive and P ridibundus often will replace 
the original population. Note that in this case we also have 
stable L-E-R complexes as the outcome. This is a system 
where two independent populations coexist, an L-E com- 
plex and a P ridibundus population. The L-E complex is 
stable due to female preferences and lethal mutations on 
the R genomes, while the P ridibundus population is sta- 
ble because of the absence of deleterious mutations, which 
are purged by selection. The two populations do not cross 
because the number of individuals in both is high enough 
to ensure that females of one population in most cases 
find a preferred male, of the same population, in the set 
randomly chosen by the best-of-n procedure. 



The whole population collapses when the introduced 
frogs have a low fitness. In this case, P ridibundus indi- 
viduals will not survive for long, given their unfitness, but 
P ridibundus frogs, before their death, can introduce R 
genomes without mutation in the P esculentus population, 
thus provoking the collapse of the complex. 

Finally, by assuming that P ridibundus are at a con- 
siderable disadvantage, the introduced unfit population is 
out-competed. During their short survival time, P ridi- 
bundus females are not able to have a sufficient number of 
matings with P esculentus males, thus they cannot intro- 
duce a sufficient number of R genomes without mutations 
into the P esculentus population. 

Conclusions 

We have presented an individual-based computational 
model to study L-E water frog populations, i.e. complexes 
composed of P lessonae and P esculentus. The individ- 
ual based model considers not only the genotypes, but 
also the age of each individual and the average lifespan. In 
addition, female preferences (implemented by a best-of-n 
procedure) and ecological selection are considered. 
We believe that our results highlight that: 

- deleterious mutations in the R genomes strongly 
contribute, together with sexual preferences, to 
maintaining the stability of such complexes, 
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Figure 9 Results of typical simulations with an initial stable L-E population in which a percentage of 5% of P. ridibundus, males and 
females, are introduced, a = 0.4, S e = 0.4, fi = 0. 
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- female sexual choice can contribute to the diffusion 
of R genomes carrying deleterious mutations in L-E 
complexes, and 

- the introduction of P. ridibundus can destroy L-E 
complexes. 
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